Magnetic Ordering and Negative Thermal Expansion in PrFeAsO 
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We report the structure and magnetism of PrOFeAs, one of the parent phases of the newly 
discovered Fe-As superconductors, as measured by neutron powder diffraction. In common with 
other REOFeAs materials, a tetragonal-orthorhombic phase transition is found on cooling below 
136 K and striped Fe magnetism with k =(1,0,1) is detected below ~ 85 K. Our magnetic order 
parameter measurements show that the ordered Fe moment along the a axis reaches a maximum 
at ~ 40 K, below which an anomalous expansion of the c axis sets in, which results in a negative 
thermal volume expansion of 0.015 % at 2 K. We propose that this effect, which is suppressed in 
superconducting samples, is driven by a delicate interplay between Fe and Pr ordered moments. 

PACS numbers: 74.25Ha;74.70.-b;75.25.+z 



Until very recently, and despite decades of research, 
no new high Tc superconductors had been found- per- 
haps due to the perceived wisdom that the only likely 
candidates were doped Mott insulators. In this con- 
text, it is therefore difficult to overestimate the im- 
portance of the discovery of superconductivity at 25 K 
in LaFeAsOi-xF^. Like the cuprates, the REFeAsO 
(RE=ra,re earth) family of compounds are layered, how- 
ever, the parent phase is a bad metal and iron is in a 
tetrahedral coordination with arsenic in contrast to the 
square planar geometry of the cuprates. On cooling, 
the undoped parent phases undergo an electronic transi- 
tion at T s ~150 K at which the resistivity and magneti- 
sation drop sharply. Initial powder diffraction studies 
showed that this transition corresponds to a structural 
change from tetragonal to orthorhombic symmetry 2 -, and 
at slightly lower temperatures, commensurate Fe-spin or- 
der develops^£. 

Initially, many authors speculated that the structural 
distortion is the result of a spin density wave transition 
and that superconductivity arises as a result of doping 
which suppresses both the structural phase transition 
and the Fe-spin ordering- 5 -. On this basis, strong analo- 
gies have been drawn with the cuprate materials. Al- 
ternatively, as the Fe-As layers in the high temperature 
tetragonal phase of these compounds are a perfect real- 
ization of the J\ — J2 frustrated square lattice, magnetic 
frustration has also been proposed as a potential driv- 
ing force for the symmetry lowering phase transition at 
150 Ki. Indeed, the magnetic structures reported for the 
majority of the undoped compounds are simple collincar 
striped models, which are predicted to be the ground 
state of the J\ — J2 model when J1/J2 ~1 [7]. Evidence 
for coupling between structural and magnetic order pa- 
rameters is particularly strong in the AFeiAsi (A = Ca, 
Sr, Ba) family of materials. Neutron diffraction studies 
show that magnetic order emerges exactly at the struc- 



tural transitions^. In the REFeAsO compounds how- 
ever, the structural and magnetic transitions are well sep- 
arated, in LaFeAsO^ and CeFeAsO 1 - 1 -, spin order is found 
some 20 K below T s suggesting a weaker coupling. 

In this communication we describe a series of neu- 
tron scattering experiments which show that in PrFeAsO 
there is a magneto-elastic coupling that results in a nega- 
tive thermal expansion (NTE) of the c— axis below ~40 K 
and an expansion of the unit cell volume. While for this 
compound T s =136 K, our measurements become sensi- 
tive to Fe magnetic order much below T s at a tentative 
TAr(Fe)=85 K, the Fe spins order in the form of antifer- 
romagnetically coupled stripes with the moments point- 
ing along the a— axis and propagation vector k =(1,0,1). 
The Pr spins are found to order below T/v(Pr)=12 K. 
The emergence of the lattice anomaly correlates with 
changes in the intensity of the (112) magnetic reflec- 
tion and is likely driven by a subtle interplay between Pr 
and Fe spin order. Our diffraction measurements show 
that the NTE behavior is suppressed in superconducting 
PrFeAsO0.85F0.15 together with magnetic order. 

We synthesised polycrystalline samples of PrFeAsO 
and PrFcAsO0.85F0.15 by the previously reported 
methods 1 - 2 -. The resistivity data for PrFeAsO shown 
in Fig. [1] exhibit several anomalies below 300 K. The 
most prominent feature is a broad maximum around 
T s . While the resistivity increases when approaching 
T s from higher temperatures there is a pronounced sup- 
pression of p around 150 K. Such a behavior is typical 
for the REFeAsO materials. Interestingly, similar to e.g. 
LaFeAsO, there is a kink in the derivative dp/dT at a 
temperature slightly below T s £ From the comparison 
with LaFeAsO we conclude that this kink indicates the 
onset of long range spin order as indirectly probed by 
the electrical resistivity, i.e T^(Fe) ~127 K. Note, that 
in contrast to LaFeAsO large RE moments are present 
in PrFeAsO which mask any anomaly in our static mag- 



FIG. 1: (color onine) (a) Temperature evolution of the 
PrOFeAs (220)< peak as measured at 2(9=82.3 deg. and 
FWHM of the (220)< and (400)/(040) o peaks, (b) Resistiv- 
ity measurements from the PrFeAsO (c) and superconducting 
PrFeAsO0.85F0.15 sample. 



netisation data around T s . An additional anomaly, i.e. 
a sharp decrease of the resistivity, is observed at much 
lower temperatures, i.e. Tjv(Pr)=llK. Our neutron data 
(see below) confirm that at this temperature magnetic 
order of the Pr spins evolves. Our fluorine doped sam- 
ple shows a superconducting transition with a Tc=43 K. 
Neutron powder diffraction (NPD) measurements were 
measured as a function of temperature on the D20 pow- 
der diffractometer— located at the Institut Lau-Langevin 
using a wavelength of A=1.88 A. For these measurements 
we collected NPD data continuously every 60 seconds 
while the sample was warmed at a rate of 2 K/min. be- 
tween 1.4 and 100 K and 1 K/min between 100 and 160 
K. Measurements on cooling between 200 and 150 K were 
made in a similar fashion with a cooling rate of 2 K/min. 
This mode of operation allows for a high degree of flexi- 
bility in the data collection and analysis. While the high 
intensity of D20 allowed us to collect high quality NPD 
data within a 1 minute period suitable for Rietveld re- 
finement, when better statistics are desired the data can 
be integrated over a desired temperature range (typically 
over a few K). Additional NPD data were also measured 
using the E9 high resolution neutron powder diffractome- 
ter with A=1.79 A and the FLEX cold triple axis spec- 
trometer with A=4.05 A, both at HZB. The NPD data 
were analysed using the Rietveld method with the pro- 
grams GSASH and FULLPROFi^. 

Turning our attention on the NPD measurements from 
D20, we find that on cooling the PrOFeAs sample, the 
structural phase transition from P4/nmm to Cmma oc- 
curs at T s =136 K, similar to what has been reported 
for the other REOFeAs material o 2 i 3 i 16 ' 17 . The evolu- 
tion with temperature of the (220)t tetragonal reflection, 
which splits on cooling, into the (400) o and (040) o re- 
flections is shown in Fig. [2ja). The sharp correlation 
between the resisitivity decrease at 136 K and the struc- 
tural phase transition is highlighted by plotting the in- 
tensity of the (220)t peak from the NPD data between 
83.3 and 83.4 deg as a function of temperature (fig. [I]). 
The exact determination of the transition temperature 
from the NPD data is challenging as the nature of this 
transition appears to be first order— <1Q. The exact tran- 
sition temperature can not be established simply by the 
sequential Rietveld refinements of the NPD data as the 
data can be almost equally well fitted just above or below 
T s with either the Cmma or P4/nmm structures. A bet- 
ter method to analyse the splitting of the (220)t reflection 
close to T s is to model it with either one or two resolution 
limited reflections. This defines the transition tempcra- 



FIG. 2: (color online)(a) Temperature evolution of the (220)< 
peak in PrFeAsO (D20 data) showing the orthorhombic split- 
ting at T s . (b) Temperature dependence of the a— and 6— axis 
(top) and c— axis (bottom) determined from Rietveld analysis 
of the neutron powder diffraction data. A clear anomalous ex- 
pansion of the c— axis, which corresponds to the layer stacking 
direction is observed below 50 K. This expansion is confirmed 
by measurement on the E9 diffractometer shown on the same 
panel as (blue) open symbols. 



ture (for example on cooling) as the point of which the 
(220)t reflection becomes significantly broader that the 
resolution of the diffractometer to support modelling by 
two peaks. Such analysis shows that the full width at 
half maximum (FWHM) of the (400)/(040) o reflection 
sharply increases at T s =136 K in excellent agreement 
with that obtained from the resisitivity measurements, 
while the (220)* reflection is signficantly broader above 
T s . A physical interpretation of such analysis would sug- 
gest an abrupt transition on warming with an extended 
range of co-existance as a broader peak within the reso- 
lution of the measurement. 

The effect of the structural distortion on the cell pa- 
rameters is also shown in Fig. (Hb). The a— axis is com- 
pressed and the 6— axis expands leading to incquivalent 
J[ioo] an d <%)io] interactions in the Fe-As planes as the 
Fe-Fe distances and Fe-As-Fe bond angles change. These 
structural changes lift the magnetic frustration, allowing 
the development of Fe spin order (vide infra) at lower 
temperature. Surpisingly on further cooling below ~ 45 
K the c— axis (which corresponds to the layer stacking 
direction) smoothly expands. As shown in Fig. [3j this ef- 
fect dominates the unit cell volume at low temperature, 
leading to a small, but measurable, negative thermal ex- 
pansion of the volume (see below) . This increase of the 
c— axis was confirmed at low temperatures also on the E9 
diffractometer at HZB in a constant temperature mode 
of operation (shown as black circles in Fig. EJc)), con- 
firming our measurements on D20. In order to estimate 
the size of the effect, we fitted the high temperature (50 
160 K) cell volume to a 2nd order Grimeiscn approxima- 
tion yielding a Debye temperature 9d = 355(5) Ki&. The 
extrapolated cell volume at low temperature is shown as 
a solid line in Fig. [3] Comparing the expected unit cell 
volume on the basis of the 2nd order Griineisen fitting 
we calculate a volume change of ca. 0.015 % at 1.4 K, 
which is remarkably similar in magnitude to the magne- 
toelastic volume effect and negative thermal expansion 
reported in other frustrated magnets, such as the spinels 
ZnCriCh\ (Ch = S, Se)^. We propose that negative 
thermal expansion in PrFeAsO originates in magnetoe- 
lastic coupling related to the onset of Fe and Pr spin 
order. Below T s we find evidence of magnetic ordering 
arising from Pr and Fe spins. On cooling below T s we 
observe the emergence of magnetic reflections, the most 
prominent being the (101) and (102). These reflections 



FIG. 3: (a) Magnetic order parameter for PrOFeAs mea- 
sured at the (102) o reflection; (b) Magnetic order parameter 
measured at the (101) reflection; (c) Refined cell volume of 
PrFeAsO from neutron powder diffraction data, red circles 
represent data collected using D20 at the ILL, blue circles are 
data collected using E9 at HZB. Also shown is a Debye fit to 
the high temperature (50 - 120 K) cell volume. 



FIG. 4: (left) Low angle region of neutron powder diffraction 
patterns for PrFeAsO, recorded using D20. Top pattern shows 
fit at 30 K in Fe spin ordered regime, bottom pattern shows 
combined Fe and Pr fit at 1.4 K. Stars mark areas of dis- 
agreement between observed and calculated patterns; (right) 
Refined Fe magnetic structure of PrFeAsO at 30 K. 



index on the chemical cell and correspond to a simple 
striped model of Fe spins with moments of 0.35(5) \ib 
along the a— axis with k =(1,0,1) at 30 K as shown in 
fig. [H At 12 K we also find that Pr-spins order giving 
a sharp rise in the intensity of these magnetic reflections 
as well as the observation of a new (100) magnetic re- 
flection. The best fit to the low temperature NPD data 
was achieved by a model where both Fe- and Pr-spins 
order within the ab— plane in the two-left, one-right mag- 
netic structure described for NdOFeAsi&iLSi. At 1.4K 
the moments of the Fe- and Pr- ions are 0.53(20) and 
0.83(9) fiB respectively. We searched extensively, but un- 
successfully for reflections which would justify a change 
in the magnetic propagation vector below 12 K, as was 
reported for RE=Ce~. In addition we attempted to fit 
the magnetic scattering at low temperature with the Fe 
spin model described above and various arrangements of 
Pr moments, with and without doubling the c— axis. The 
model described above provided the best fit to the data 
by far. 

The onset of the Pr and Fe spin ordering was deter- 
mined by following the temperature dependence of the 
magnetic reflections using both the D20 data as well as 
performing additional measurements on the FLEX spec- 
trometer. The temperature dependence of the (101) and 
(102) magnetic reflections are shown in fig. ^ a) and (b). 
Given the small size of the ordered Fe moment, the D20 
data were binned in approximately 5 K intervals to en- 
hance the statistics of a given diffraction pattern. To ob- 
tain the temperature dependence of the (102) reflection 
the data were integrated over a 1.5 deg range around the 
peak. A similar size range of 28 adjacent to the (102) 
reflection was also integrated to gain an estimate of the 
background and its temperature dependence, which was 
found to be constant. The results of this analysis are 
plotted in fig. 02a). We find that although there is a well 
defined Bragg reflection at low temperature the mag- 
netic scattering becomes much weaker as temperature 
is increases towards T s . A sharp drop in the intensity 
of this reflection is observed on warming through 12 K, 
which signifies the onset of Pr ordering, and the magnetic 
scattering finally disappears above 85 K. This observa- 
tion will tentatively place Tjv(Fe)~85 K, however, In the 
present investigation we are only sensitive to magnetic 
scattering once the Fe moment reaches a value of 0.1 he- 
more sensitive probes to Fe-magnetism such as muon and 
Mossbauer spectroscopy may well place Tjv(Fe) to be ei- 



ther coincident with T s or are a few Kelvin below as it 
has been found for RE=L&, Ce and NdSJiSi. 

We followed the temperature dependence of the (101) 
reflection with the FLEX triple axis spectrometer and 
found its evolution to be somewhat different from that of 
the (102). A large decrease is also observed for this re- 
flection on warming through Tjv(Pr) , however, the scat- 
tering from the Fe spin order reaches a plateau around 
40 K before disappearing above 75 K. The apparent 
lower Tjv(Fe) given by this reflection likely arises from 
its smaller magnetic structure factor compared to the 
(102). 

As is evident in Fig. [3] the expansion of the unit cell 
volume coincides with the plateau of the (101) reflection 
suggesting that its origin is magnetoelastic. The rise in 
intensity of the (101) Fe magnetic reflection up to ^40 
K and its decrease as T w (Pr) is approached suggest a 
possible interplay between Pr and Fe. One possible ex- 
planation is a decrease in the ordered Fe moment along 
the a axis due to the influence of the Pr moments. Such 
behavior has also been alluded to in ref.— for CeFcAsO, 
however in the present case the magneto-elastic coupling 
that is generated is sufficiently strong to produce an in- 
crease of the unit cell volume. Such effects have not been 
observed to date in the REFeAsO family of compounds. 

We also performed neutron powder diffraction mea- 
surements on superconducting PrFeAsO0.85F0.15 down to 
0.4 K using a 3 He insert. We detect no evidence for Fe 
or Pr magnetic order, while the NPD data show a posi- 
tive thermal expansion for the a— and the c— axis. This 
result is also in contrast to measurements for cuprate 
superconductors such as Nd2- K Ce a; Cu04, in which rare 
earth antiferromagnetism is found to co-exist with the 
superconducting stated below 2.5 K. This would indi- 
cate that the rare earth magnetism in the layered FeAs 
compounds may be induced by the Fe-spin order. 

In conclusion, our measurements show that PrFeAsO 
undergoes the same T - O phase transition as the other 
REFeAsO compounds at 136 K. On cooling below ^85 
K we detect long range magnetic order of the Fe mo- 
ments in a striped arrangement, with k =(1,0,1) and at 
lower temperatures (12 K) long range order of Pr spins. 
Our detailed Rietveld refinements have identified a previ- 
ously unknown coupling of magnetic and lattice degrees 
of freedom at ~40 K, which leads to an anomalous ex- 
pansion of the layer stacking direction. We propose that 
that this negative thermal expansion results from a deli- 
cate interplay between the ordered Fe and Pr sublattices. 



Suppression of the NTE state and long range magnetic 
order is necessary for the emergence of superconductivity 
in Pn-^FajFeAsO. 
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